FERMI SPECTRA AND THEIR GAUGE INVARIANCE IN HOT AND 
DENSE ABELIAN AND NON-ABELIAN THEORIES 
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Abstract 



The one-loop Fermi spectra (one-particle and collective ones) are found for all momenta in 
^ ; the T 2 -approximation and their gauge invariance in hot and dense Abelian and non-Abelian 

theories is studied. It is shown that the one-particle spectrum, if the calculation accuracy is 
^ ■ kept strictly, is gauge invariant for all momenta and has two branches as the bare one. The 

collective spectrum always has four branches which are gauge dependent including also their 
|q| = limit. The exception is the case m, // = for which this spectrum is gauge invariant 
<3\ ' for all momenta as well. 
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1 Introduction 



The problem to obtain the gauge invariant results from gauge theories (especially from non- 
Abelian ones) is not new but is very actual for present-day physics. This problem is very 
many-sided and one easily calls many tasks which require the first rate solution. In particular, 
many questions arise when the effective electron (or quark) mass is calculated within hot and 
dense QED or QCD, including even the pure QFT (quantum field theory). Of course, the 
latter case is more studied but some problems (in the first turn, the proof that the fermion 
mass is infrared-finite) were solved only recently making calculation of the perturbative 
mass in QFT selfconsistent and proving its gauge invariance [1]. The completely different 
situation takes place when T (temperature) is nonzero. In this case the nonzero thermal mass 
(or thermal gap) is always generated [2,3]) and all Bose and Fermi spectra (if we speak about 
so called collective ones) are split demonstrating themselves more pronounced than in QFT 
[2-9]. The more complicated scenario arises when the bare fermion mass m and temperature 
are nonzero simultaneously [10,11]. In this case the spectrum has the additional splitting 
and changes its long wavelength asymptotical behaviour: it becomes q 2 instead of |q| when 
m = 0. If along with m,T^0 the chemical potential \x is also nonzero the scenario becomes 
very cumbrous [12,13] and only the separate spectrum limits can be found analytically. 
To obtain the spectrum curves for all |q| the numerical calculations are necessary that, 
unfortunately, hides many details from the further analysis. However, if scenarios with 
/i, T 7^ and with m, T ^ are considered separately the spectrum curves for all |q| 
can be found analytically in any case within T 2 -approximation [14,15]. These analytical 
expressions easily reproduce all known limits and are more convenient to investigate the 
spectrum gauge invariance and many other properties. For the general case when /i, m along 
with T are nonzero only the effective thermal mass can be found analytically [15,16] and 
being calculated within different gauges at once demonstrates its gauge dependence. However 
the situation drastically changes if the one-particle spectrum [16] is considered instead of the 
collective one. This spectrum is analogous to the bare one but its properties are modified 
due to interaction with medium and completely different from the collective ones. Unlike the 
latter scenario the one-particle effective mass is gauge invariant in the leading e 2 -order (and, 
possibly, in the higher orders as well) and has a real physical sense. The nonperturbative 
fermion mass (the |q| = 0-limit of the collective spectrum) is always gauge dependent as well 
as the full collective spectrum, and the additional resummation seems to be necessary. 

Briefly speaking the question of the gauge invariance for many results found in statistical 
QCD is insufficiently studied and opens for discussion. Even in statistical QED the exact 
gauge invariance is only known for the photon thermal mass, but the same as in QCD, no 
exact results exist for the thermal fermion masses and each times their gauge invariance is 
necessary to investigate independently. 

The goal of the present paper is to show that not only the perturbative mass but also the 
full one-particle Fermi spectrum for all q 2 is gauge invariant if the calculation accuracy is 
kept strictly. This invariance takes place for any m, /x-parameters, although the case /i = 
is separated. In the case /i = the full one-particle Fermi spectrum is gauge invariant even 
beyond the perturbative accuracy in any case within the simplest summation of one-loop 
calculations. In the section 4 the collective spectrum is calculated for the case m, fi ^ in 
the Coulomb gauge to compare it with one found in the Feynman gauge. It is shown that this 
spectrum is always gauge dependent but, the Coulomb gauge seems to be more physical one 
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for any applications. The exception is the case to, \i = [3,5] for which this spectrum is gauge 
invariant for all momenta including their effective mass. This scenario, probably, is valid in 
all perturbative orders since in this case within T 2 -approximation no dimensional parameters 
(except T) are present and due to general theorems [17,18], such spectrum should be gauge 
invariant exactly (including all higher order corrections) even within hot QCD. Of course, 
the Abelian QED and non-Abelian QCD should generate the different gauge dependence of 
any spectra, but in the leading e 2 (or g 2 ) order QED and QCD Fermi spectra are the same 
(if everywhere the numerical factor is changed due to the prescription e 2 — > g 2 (N 2 — 1)/2JV). 
Today there is a problem to find (or, at least, understand) what a difference arises within 
the spectrum found in QED when the higher order corrections are taken into account to 
distinguish QED from QCD. The strong infrared divergencies of hot QCD (see, e.g. [19] 
and other references within it) should display themselves via the higher order corrections 
and it is not excluded that effective thermal mass in QCD becomes the infrared-unstable. 
To start we choose hot and dense QED (not QCD) although all results are the same in the 
approximation considered. We also propose that everywhere the damping is small and can 
be neglected considering that this question should be investigated separately [20,21]. Our 
attention will be focused to solve the fermion dispersion relation in a more complete form 
and to investigate the gauge invariance of all spectra found. 

2 The one-loop fermion self-energy in the Coulomb 
gauge 

The Coulomb gauge is rather reliable gauge for perturbative calculations within hot gauge 
Abelian and non-Abelian statistical theories. It does not generate (unlike the a-gauges) the 
additional infrared divergencies and requires only the standard ultraviolet regularization. Of 
course, all calculations in the Coulomb gauge are more complicated than in the Feynman 
one but these are the technical difficulties which can be overcome without any principle 
modifications of the theory studied. In what follows our calculations are performed in the 
Coulomb gauge but often we compare them with the ones made in the Feynman gauge [14- 
16] and keep the same abbreviations. Some details of these calculations and their cumbrous 
results are placed in Appendix A. 

Within the one-loop approximation the exact decomposition for S(g) is given by 

E(q)=i lll K li (q)+mZ(q) (1) 

where three new scalar functions are introduced to find nonperturbatively the function G(q) 

r( v = -h^ + K^) +m (1 + Z) 

{q) (q, + K,y + m 2 (i + zy ■ [) 

This representation leads to the one-loop dispersion relation (the exact one has the additional 
functions [22]) for the Fermi excitations which in any gauge has the form 

[ (iq 4 - /i) - K A } 2 = q 2 (1 + K) 2 + m 2 (l + Z) 2 (3) 

and after the standard analytical continuation it can be solved analytically or numerically. 
Here K± = %K± and q — {(g 4 + i/i), q}. 
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The calculations are rather lengthy and require at first to extract the functions Z(q) and 
K^q) from 

\K 74£ P + ij(p - q) [(p - q|p)/(p - q) 2 ] + m 



s c (g) 



f d 3 P f r< 

7 (2vr)3 \ U p 



[g4 + «(/i + e p )] 2 + (p-q) 2 
n|_ (|p| +fi- iqihi - ?7P [(p|p + q)/p 2 ] - "i 
Pi 



[g 4 + i(/i+ |p|) P + Cp+q 

™p 74e P - hp - m 
e 



h.c; (m, /i) — > — (m,/i) 



(4) 



p 2(p-q)2 

and then to calculate them using a number of exact integrals. Our result for the functions 
Z(q) and K 4 {q) has the form 

oo 



Z c (q 4 , q ) = e 2 J 



d\p\ |p| j 1 \n+ + n p 



4tt 2 2|q| ( e p 



( ln(a£a F ) 



-In 



(|p| + |q|) s 



iFT 4 c (g 4 ,q) =e 2 

(|p| + |q|) 2 
(|p|-|q|) 2 



- |q|) 2 

d|p| |p 



ln(a£a B ) 
P 



(5) 



In 



4tt 2 2|q| 



ln(^' 



(fi - iq 4 ) 



( ln(a£a F ) 



]n(a£a B ) 



(6) 



and the more complicated calculations are necessary to obtain the vector K n (q) where n = 
1,2,3. Here we use a definition K n (q) = q n K(q) and after that the scalar function K(q) is 
calculated to be 



q 2 ^ C (g4,q) = e 2 
+d F \n(apap)) + 



djp[ f P 2 
2tt 2 I e n 



n+-n p 



1 



n+ + n p 



1 + 



8|p||q| 

1 a + 
- — n — r( hp \n(ata F ) + dp ln(— 

|p||q| v a F 



1 a + 1 

+nj |p| [1 + ——(h B \n(a+a B )+d B ln(^) ) ] + q 2 X c ( ?4 ,q 
o|p||q| a B J 



) 



(7) 



where hp = q 2 — "m 2 — (i?4 — At) 2 and d F = 2e p (ig 4 — /i) (analogously hs = q 2 + m 2 — (ig 4 — 
yu) 2 and ds = 2|p|(ig 4 — //)) are some simple abbreviations and other ones have the more 
complicated form 



q 2 — m 2 


- (m 


- fi) 2 ±2e p (iq A 


-^)-2|p||q| 


q 2 — m 2 


- (m 


- fi) 2 ± 2e p (ig 4 


-A*) + 2|p||q| 


q 2 + m 2 


- (m 


-/i) 2 ±2|p|(^ 4 


-^) - 2 lp||q| 


q 2 + m 2 


- (m 


-/x) 2 ±2|p|(*g 4 


-/i) + 2|p||q| 



(8) 



(9) 



Here and everywhere e p = v^ 2 + P 2 - The function K c (q4, q) has the rather cumbrous form 
and is placed in Appendix A. This function is mainly essential beyond the T 2 -approximation 
and below only a few of its terms are exploited to find /i/T-corrections. 
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3 The one-particle spectrum of massive Dirac particles 
in the T 2 -approximation 



The one-particle spectrum is the perturbative one and corresponds to the calculations when 
in the leading order the dispersion relation (3) is solved with £(§4, q) taken at once on the 
bare mass shell iq 4 = // ± y/q 2 + m 2 . 

Within e 2 -approximation this spectrum can be presented as follows: 



iq A = (» + K 4 ) ± yj w? (l + Z{q)) + q 2 (l + K{q)) (10) 

where all functions being put at once on the bare mass shell are independent on iq A . This 
spectrum is qualitatively different from the collective one and at the beginning is more useful 
for many applications. In particular, its long wavelength limit (|q| = 0-limit) which defines 
the effective thermal mass is the gauge invariant value [16] 

iq 4 = n R ± m R = /I ( 1 + 2 I B ) ± f m(l - 41 z ) + ^ 1 (11) 

At any rate Eq.(ll) is the same in both gauges: the Coulomb and Feynman ones. Here 
I B = jiI B and other abbreviations are: 



00 

j - T F, T B_ p 2[ d\P\ n p +n p r d|p| B 

Ik - I K + Ik - e J ^ e p + e J \p\n p 



00 - 1 I I _ 00 -I I I _ 

t _ 2 fd\p\ 71+ -n p _ 2 rd\p\ 7i + + 7i p 

p 

The found gauge invariance is a very important property of this spectrum and our task is to 
investigate whether this property keeps for all momenta or only some separated spectrum 
limits demonstrate it within the T 2 -approximation. However, checking this property one 
should be careful exploiting Eq.(10), since it is given beyond the accepted one- loop accuracy. 
The more reliable expression has the form 

«, 1 = (,± f ,) + A' l(q )±(!^ + 3^M) (13) 

^q 

where the one-loop corrections to the bare spectrum are given strictly within the e 2 -accuracy. 
Below we return to this question once more to demonstrate the important result: the gauge 
invariance is a priori broken if the perturbative series are considered beyond the accepted 
accuracy. 



3.1 The fermion one-particle spectrum with m,fi ^ in Coulomb 
and Feynman gauges and its gauge invariance 

Here the one-particle Fermi spectrum will be found for all momenta in the most general case 
m,/i^0 using the expression for E(q) obtained in the Feynman (F.G) and Coulomb (C.G) 
gauges. Namely these gauges are preferable in statistics since only they are reliable without 
any additional regularization. 
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Going to the T 2 -approximation, we simplify all logarithms which define the above found 
integrals, keeping only the leading T 2 -term and /i/T-corrections. Effectively this operation 
leads to the following ansatz 



ln(-£) ~ £*(q) 



1 



2q 



rn 



1 + 



2q 2 



;i±fi) 

q 



ln(=#) ~ £f (q) = In 



|q| 
|q| 



(14) 



and also ln(a£a F ) ~ — 4|q|/|p| and ln(a^a B ) ~ 0. In Eq. (14) and everywhere the bottom 
signs of £±^ B -quantities correspond to the signs in the bare spectrum iq± — /j, ± y/q 2 + m 2 
which is inserted into all integrals before any expansion. The derived ansatz allows to 
calculate the T 2 -terms and /x/T-corrections for all functions involved in Eq.(10) explicitly 
and to compare the spectra found in different gauges. The results of the calculations made 
in the Coulomb gauge 



Z c (q) = -U z - 



2|q| 



^f(q) 



4|q| 



f ^(q) B 

w K 



T r - 



q 2 ^(q)=/^ Teq ir 4 C '(q) 
and in the Feynman one 

Z*(q) = Z*(q) - ^M/ B 
v 4|q| 



(15) 



Ki"(q) = iff - 2/b 



q 2 iTf (q) = 2ml B 



±^) + q 2 ^(q) 



m 



(16) 



are very similar but are not equal if // ^ 0. This means that £((74, q) and its structure 
functions rule, gauge dependent but this is not the case for the one-particle spectrum 

since within Eq.(13) the algebraic cancellations are possible to restore the spectrum gauge 
invariance. 

Indeed within the Feynman gauge the one-particle spectrum, before algebraic cancella- 
tions, has the form 



iq A 



(/x ±Cq ) + (*f(q)-2J fl ) ± {^(Z c (q)-^/ 



+ — 



2m ^(Trf ±-) + q 2 ^) 

v 4 q m 7 



4|q 



(17) 



but can be easily simplified as follows 



(/i±6 q ) + K 4 C (q)± ( 



m 2 Z c (q) q 2 ^(q 



+ 



(18) 



that is exactly the one-particle spectrum in the Coulomb gauge, and consequently the gauge 
invariance is restored. 
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However the cancellations are not complete if the one-particle spectrum is presented by 
Eq.(lO). In this case the spectrum found in the Coulomb gauge 



<?4 = (jm + ^ 4 c (q)) ± y™ 2 (i + ^ C (q)) 2 + q 2 [i + ^(Mq) =f (q))] 2 (19) 

is essentially different from the one found in the Feynman gauge 
^ 4 = (/i-2/ B + i? 4 c (q))± (20) 



m\\ + ZO(q) - ^/b) 2 + q 2 [l + ^(^(q) =F (q) + 2mI B (^M ± fa)) 
\ 2|q| L q 2 v 4|q| m ' 

since the accepted accuracy is exceeded within this scenario. These spectra are gauge de- 
pendent although their low and high energy limits reproduce the gauge invariant results as 
well. The correct one-particle spectrum should be found using Eq.(13) and has the form 

iq* = (M±*d±[ - ^(4/ z+ ^Ml B ) + **] (21) 

e q ^|q| e q 

which is our main result for this section. To find its limit for small q 2 one should use the 
standard expansion of £± B (q) -quantities 

that is the same for both functions in the leading order and perform the simple algebra. 
These calculations yield the following result: 

*g 4 = /i(l + 2/ fl ) ±m F ± + + 0(q 3 ) (23) 

where trip = m(l — AIz) + Ik /mis the effective fermion mass accepted for e 2 -approximation. 
The high energy limit in the T 2 -approximation is also easily calculated and has the form 

iq A = pi ± |q| ± l-[I K -Am 2 I z ] + 0{\\n^-) (24) 
|q| q 2 2m 2 

which is analogous to the bare one. All other terms are small and can be neglected. Com- 
paring the low and high energy limits we see that one spectrum branch (the case of an upper 
sign ) can demonstrate the minimum at finite momentum if the density is rather large. 

3.2 The gauge invariant one-particle spectrum with fi = in the 
T 2 -approximation 

The special case of zero density (when /i — 0) is rather interesting since in this case within 
the T 2 - approximation there is no problem with the gauge invariance at all. In this case the 
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one-particle Fermi spectrum is gauge invariant at the beginning and Eq.(lO) can be used to 
present it as follows 



iq 4 = Kt( q ) ± J m 2 (l - 4/z) 2 + q 2 [ 1 + ± ( I K =F ^(q) ) ] 2 (25) 



where 



«f (q) = -f¥ '£ -If'^^Tiri + °(i 4 ' < 26 > 

4|q| 4|q| m 6m 6 

Its low and high energy limits repeat the appropriate expressions found above in the Coulomb 
gauge for the case \x = but, and it is more important, this spectrum is gauge invariant at 
once and, probably, this property is kept within all higher order corrections 

4 The collective Fermi spectrum and its gauge depen- 
dence in the T 2 -approximation 

For hot gauge theory (Abelian or non-Abelian) with massless fermions and in the symmetrical 
case (/i = 0) the Fermi excitations have collective nature and appear as "massive" quasi- 
particles (or quasi-holes) under the physical vacuum. In the simplest case their thermal 
masses (or more exactly their thermal gaps) are the same for each spectrum branch and 
they arise dynamically in spite of the exact chiral symmetry inherent to initial Lagrangian 
on the operator level. Such excitations in hot QCD were first found in [3] and then studied 
in many other papers. Their spectrum [3,5] was found analytically in the T 2 -approximation 
for all momenta 

4(0 = e 2 - 2 (^ + |ln|^), V = ±l (27) 

and is gauge invariant in any case within one-loop calculations. Here the variable £ runs over 
the range 1 < £ < oo and the long wavelength spectrum limit (i.e. q = 0-limit) corresponds 
to the case £ — > oo. In QCD the thermal gap is: Uq = g 2 T 2 /6, and the spectrum is split since 
rj = ±1. All its branches have the finite gap at zero momentum (the finite thermal mass) 
and a> 2 (£) quasi- hole branches have a very specific minimum at finite q. This minimum 
always exists when m, \x = and is a very interesting subject for discussion [3-9]. For small 
q this spectrum has the linear limit 



2 |q| 7 q 



2 



u; 2 (q) = <4 1 ± - if + - ^ + 0(|q| (28) 

since all branches have the finite gap, but no real mass. In the high momentum region the 
quasi-particle and quasi-hole branches have qualitatively different limits. The quasi-particle 
branches demonstrate the standard powerful behaviour 

c4(q) = q 2 + 2ul - (^) Ini + O(^) (29) 

q (jJq q 
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but this is not the case for the quasi-hole ones 




which more quickly (exponentially) approach the line uo = q. 

However this scenario changes drastically when m or \x is not equal to zero. In the first 
turn this spectrum becomes gauge dependent and more complicated: its mass (or gap) is 
split and its minimum inherent to the quasi-hole branches, as a rule, disappears. We study 
this situation below using the Coulomb gauge and compare our calculations with the results 
found in the Feynman one from the papers [14,16] where these excitations were also obtained 
analytically for all momenta with nonzero m, /i-parameters. 

4.1 The collective spectrum of massless Dirac particles in the 
Coulomb gauge for the case fi ^ 

The scenario with m = but /i ^ is very similar to the previous one where m, \x = and 
keeps many its properties although the symmetry {iq± — > — iq±) is lost. Here, as previously, 
all spectrum branches have the finite gap and the linear limit for small q, but there are two 
different gaps one of which, sometimes, can be equal to zero. 
The dispersion relation (3) is simplified to be 

[ (jg 4 _ ^ _ tf 4 ]2 = q 2 + (31) 

and its solution can be written as follows 

(iq 4 -p)-K 4 = r]\q\ (1 + K) (32) 

Here rj = ±1 and we use the simple redefinition K A = iK±. All functions involved into 
Eq.(32) are calculated using the T 2 -approximation where the essential simplifications are 
possible due to the simple ansatz 

Ma + a~) = -HM + O(-L) , In £ = 2 In + O(-L) (33) 

which keeps only T 2 -terms and /i/T-corrections in Eq.(32). With the ansatz (33) the further 
calculations are easily performed and their results are given by 

^ 2 V ^ J f + 1 / |q| 

(34) 

where all numerical integrals Ik and are defined by Eq.(12). Now one should plug the 
expressions found above into Eq.(32) and perform a simple algebra to find uo = £|q|. Here 
uj = {iq± — /i) and the variable £ is more convenient than |q|. The latter should be excluded 
from Eq.(32). The result is the quadratic equation with respect to cj(£) 

J- (£-,) - w(Ib( 1-1)6(0) " I K eA(O=0 (35) 



^ (94 , q) = ^( 1 + ! ln i^)- 

i?f( 9l ,q) = / B -^lnf^ 
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where the functions A(£) and &(£) are given by 

KO = e-^(l-ain|^ (36) 
Keeping the e 2 -accuracy our solution of Eq.(35) is found to be 

= ? I 2({-,) V f--^ + 2 ^~+T ' I (37) 

and presents the one-loop spectrum in medium for the case /i ^ 0. Here 77 = ±1 as usual 
and uj = iq 4 — /i. The spectrum (37) has four branches which are split for all momenta, 
excepting their |q| = limit, where they are combined by pairs demonstrating two different 
gaps 

E ± = ^ ± sJT K (38) 
Their asymptotic behaviour for small |q| has the form 

u ± (q\ V ) = E ± + I |q| + \ ( I B ± I JT K ) i + 0(q 3 ) (39) 

where one branch always has the minimum at finite momentum. Two branches, if E_ = 0, 
lose their gap, but this can occur only under the special constrain on /i, T-parameters. 

Unfortunately this scenario is gauge dependent and spectrum (37) differs from the one 
found in the Feynman gauge [14]. Nevertheless, some mathematical correspondence takes 
place: the spectrum in the Feynman gauge arises from Eq.(37) if u;(£) will be replaced to 
[— u;(£)]. This correspondence is rather strange and is inherent only to this simple scenario 
where the mixing is absent between \x and m parameters. If //, m 7^ the situation is more 
complicated and the gauge invariance in this case seems to be very problematical. 

4.2 The collective spectrum of massive Dirac particles for the case 
fi = in the Coulomb gauge 

It is another scenario which can be considered analytically for all |q| in the T 2 -approximation. 
Here Ib = identically and symmetry (ig 4 — > —iq^) is restored. But now the dispersion 
relation has the additional mass term 

[ {iq 4 - /i) - K A } 2 = q 2 (1 + Kf + m 2 (l + Z) 2 (40) 

and all calculations are more complicated. To solve Eq.(40), taking into account only the 
leading m/T-corrections, we use the functions from Eq.(34) with /j, — 

R°M = §(i + |hf^r). ^•i) = -^f 1 "|TT (41) 
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and Z G (g 4 ,q) = — 3I Z . With these functions Eq.(40) is transformed to be 



o/(£ 2 -l) - u 2 e[m 2 (l-3I z y + 2I K ] - Ji£ 4 (i + £hi 



= 



(42) 



and becomes the quadratic equation in respect to u; 2 (£) (not to u;(£) as previously). However, 
this is not the case when m,/j ^ simultaneously (see [15] and Appendix B). In this case 
the dispersion relation generates the full equation of the fourth degree in respect to u;(£) and 
requires the numerical calculations. The analytical solution found within Eq.(42) has the 
form 



MO 2 = 



K 



m 



R 



2(e - 1) 



± 



m)lK) 2 + m R (I K + m R /±) 



(43) 



and presents the collective Fermi spectrum for all |q| in the Coulomb gauge. Here m R = 
m 2 (l — 3Iz) 2 and, namely, this quantity differs the Coulomb and Feynman gauges, where 
m 2 R = m 2 (l — 2I Z ) 2 - Excepting this quantity the found spectra are completely the same in the 
both gauges, and, in some sense, are gauge invariant. It is not excluded that the additional 
resummation is necessary to improve the situation, but today this ansatz is unknown. The 
found spectrum is rather complicated and presents four completely separated branches: two 
branches (when the plus sign is taken in Eq.(43)) correspond to quasi-particle excitations and 
the other two (when the minus sign is taken) to quasi-hole ones. These spectrum branches 
differ in their asymptotic behavior and in many other properties. 

The long-wavelength behavior of these spectrum branches (when £ — > oo) has the form 



^(|q|) 2 



M 2 ± 



+ (m 2 



± 



I 2 



K 



m R [m R 



+ 4/ K ) 



Ml 



+ 0(|q| 4 ) 



(44) 



where the squares of the effective masses are given by 



Ml 



m 



R 



+ Ik± 



rn 



2 ( m R 



+ Ik) ■ 



(45) 



These masses are different for all branches M± = ^(r]m R ±ym R + 4/^-) and are in agreement 
with the results found in many other papers. Here rj = ±1. It is also important that the 
quasi-hole branches c<j_(|q|) 2 are very sensitive to the choice of the parameters m, T. In many 
cases these branches are monotonic functions for small |q| 2 , and the well-known minimum 
[3] disappears. Although this minimum always exists for massless particles, the special 
conditions are necessary to generate it when m^fl. 

In the high-momentum region the asymptotical behaviors found for the quasi-particles 
and quasi-hole excitations are also completely different. The quasi-particle branches are 
approximated within a rather usual expression 



«M|q| 



|q| 2 + (2I K + m 2 R ) - 



In 



4|q| 



21 



K 



m 



R 



1 



+ o(- 2 ) 



(46) 



where the nonanalytic term is not essential. The situation is different for the quasi-hole 
excitations, which do not exist in the vacuum (when T and ji are equal to zero). They 
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disappear very rapidly, and their asymptotic behaviour is found to be 

/: ix2 I 1 2 r -i , / |q| 2 (2/^ + m%). / _. . 2\q\ 2 (2I K + m 2 R ). . /{ _. 
^-(|q|) 2 = |q| 2 1 + 4 exp(- IHI *L) + 0( exp(- lHl *>-) ] (47) 

1 K 1 K 

In the high momentum region these spectrum branches approach the line oo 2 = |q| 2 more 
quickly than (46). 



4.3 The collective thermal mass in the Coulomb and Feynman 
gauges and its gauge dependence 

In the most general case when /x, T and to are nonzero simultaneously, the effective mass 
can be also calculated analytically [14-16]. Here we compare it with the results found above 
and discuss its gauge dependence. In the Coulomb gauge this mass (see [16] and Appendix 
B) is given by 



uZ(p) = -[r t m% + I B 



± 



[77 mg + I B 



+ (I K + 2rjmI B ) . 



(48) 



where m R = to(1 — SIz)- The mass found is split and has the well-separated mass spectrum, 
where two values of them pertain to the quasi-particle excitations and other two present the 
quasi-hole ones. This is evident from Eq.(48) where 77 = d=l. However, this spectrum is not 
the same as in the Feynman gauge [14] 



<4(0) = 2 



77 m F R 



± 



[77 m£ - I B 



+ (I K + ArjmI B ) ■ 



(49) 



where = m(l — 2Iz) and the /i-dependence is different. Thus, these masses, excepting the 
case to, fj, — 0, are gauge dependent and this is true although Eqs.(48),(49) are given beyond 
the accepted calculation accuracy. If one strictly keeps the e 2 -accuracy these equations are: 

u c ± = ifa ± l)m c R ± ^ ± ^(577 ± 1) , u F ± = 1(77 ± l)m F R ± ^ ± ^(7^ 1) (50) 

where we propose that to 7^ 0. If at once to = one has the simpler result 

f f - Ib 

which repeats Eq.(38). Eq.(50) leads to the following mass spectrum 

Ib 



± 



2 



(51) 



TO 



(C/F) 



V ( ™ R 



(C/F) 



- —) + M B 

m 



(52) 



m c h = 77 ^ - 21 B , m F h = 77 ^ - 41 B (53) 

TO TO 

which demonstrates that the gauge invariance for the quasi-hole thermal masses in (53) is 
destroyed more sharply than for the quasi-particle masses in (52). It is very likely, that in this 
situation the futher resummation easily restores the gauge invariance of the quasi-particle 
thermal mass even if \x 7^ but the problem with the quasi-holes one, probably, remains. 
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5 Conclusion 



To summarize we have established two kinds of fermion spectra in statistical QED and the 
same in QCD: the one-particle spectrum and the collective one. The one-particle spectrum 
is qualitatively the same as the bare one and is gauge invariant, in any case within the 
one-loop approximation. It has two branches iq^ = /ir ± jjir and their chemical potential 
and effective mass are only quantitatively changed due to the interaction with the medium. 
It is not excluded that within QED this spectrum is gauge invariant in all higher orders as 
well, if the calculation accuracy is kept strictly. Moreover in hot QED this spectrum seems 
to be gauge invariant exactly but this is not the case, a priori, for hot QCD where the strong 
infrared divergencies spoil this scenario. 

The collective fermion spectrum was calculated in the Coulomb gauge and was compared 
with the one found in the Feynman gauge to investigate its gauge dependence. This spectrum 
is additionally split (usually for all momenta) and its branches always have the nozero thermal 
masses different for all the spectrum ones, if m ^ 0. These masses arise dynamically and 
are nozero even if the exact chiral symmetry forbids them on the operator level. They are 
not correlated with bare mass and are generated always: but only for the case m, \i = this 
thermal mass and the full collective spectrum are gauge invariant. For any other cases the 
collective spectrum is gauge dependent and its connection with the real excitations should 
be investigated separately. 

However one can see (comparing our results in the Feynman and Coulomb gauges) that 
gauge invariance is broken variously in different scenarios. This fact demonstrates itself more 
sharply when /i ^ [due to asymmetry (ig 4 — > — iq^)} and is smoothly expressed (only slightly 
shifts all spectrum branches) for the symmetric case where \i = and uj 1 (not uj) presents 
all the spectrum ones. In the last case (namely, when /i — 0), the additional resummation, 
probably, can change the situation to restore the gauge invariance. But this is very unlikely 
for the scenario with [i ^ since in this case the spectra found (especially the quasi-hole ones) 
are rather different and their rearrangement seems to be difficult. Moreover this problem is 
more aggravated if all calculations are performed beyond the T 2 -approximations. In this case 
one at once encounters many difficulties arising due to the ill-infrared behaviour inherent to 
hot QCD (or QED as well) and all calculations become very ambiguous. The coexistence of 
the infrared peculiarities of hot gauge theories (like QCD) with the gauge invariance of their 
results is very problematical and this remains a very serious problem for any calculations 
performed today (especially for nonperturbative ones). 

The additional problems arise if the damping is taken into account [20,21], especially 
when these calculations are performed beyond the T 2 -approximation. Here the selfconsistent 
calculations are only acceptable, where, in the first turn, the calculation accuracy is strictly 
kept, and the resummation scheme is checked to be, at least, gauge covariant (the Ward 
identities should be satisfied between the effective propagators and vertices). 
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APPENDIX A 



The one-loop fermion self-energy in the Coulomb gauge, after the standard summation over 
the spinor indices and frequencies is presented by Eq.(4). Its calculation strongly differs from 
the calculations performed in the Feynman gauge only on the last stage when integration 
over angles is performed. This difference demonstrates itself more pronounced only for the 
vector K n (q 4 , q) which has the form 



2e P (p - q) 2 



(54) 



+ Pn ( 1 + 



pq \ n p 



p z ' ipI V («p + 



) 2 + 6 p+q 



n 



+ 



+ 



-ip + q 4 ) 2 + e 



2 

p+q 



+ 



(p-q) 27 e p V (^ P + g 4 ) 2 + (p-q) 2 ( -»e p + q 4 ) 2 + (p - q) 2 J J 

and we propose that lf„(g4,q) = (foX (q4, q). Then the very lengthy but standard calcu- 
lations lead to the result presented in Eq.(7) where the function K°\q 4 , q) is given by 



q 2 K c (q 4 ,< l ) = I 1 + I 2 + I 3 , 



(55) 



/ 



d\p\ P 2 ^ 



1 



4tt 2 



e p 2 4| P ||q|( ep + g 2 ) 2 
(P + q) 2 



-(p 2 - q 2 ) 2 dp In 



(p-q) 



, + hp dp (-e 2 + q\ ) ln(a£a F ) - dp In 

G F 



+ g( ^ + 4-4p'q 



e 2 + g 2 ) In ^ - d F ln(a^a F ) 
a F 



h = -e 1 



d\p\ p 2 ™+ + nl 



4vr 2 



2 4| P ||q|( e 2 +g 2 ) 2 



{[(-4 + ^)(p 2 -q 2 ) 2 



+ (4 + g 4 2 ) 2 (p 2 + q 2 ) 



fp — h C|) t T / 2 ^2 \ i 7 / 4- \ 

In- — + ftjr d F (-e D + g 4 ) m— - d F \n(a F 'a F ) 

(p-q) 2 a F 



+ -( /4 + 4 - 4 P q ) (-4 + ?4 ) ln(a>^) - dp In 
2 a p 



h = 



d\p\ < 

47r 2 |p| 



h B + 



P 



(h 2 B + d 2 B ) \n(a B a B ) + 2d B h B In 



JL 



The abbreviations are the same as in Eq.(7). The integrals found are essential for the next- 
to-leading orders beyond T 2 -approximation. Here only the first integral in Eq.(55) is useful 
to find /i/T-corrections. 
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APPENDIX B 



Our starting point is the dispersion relation (3) 

[ (iq 4 - fx) - K A } 2 = q 2 (1 + Kf + m 2 (l + Z) 2 (56) 

with m 7^ and we solve it in the T 2 -approximation with ix ^ 0. All calculations are 
performed in the Coulomb gauge taking into account only the leading T 2 -terms and fi/T- 
corrections. In this case the functions within Eq.(56) are: 



- K A (q A , q) 



— — r In — 1 B 

2|q| C + l 



-Zfo.q) = 3/ z + ^ln|— 4 (57) 



making it possible to essential simplify its solution. Here the variable £ = oo/\q\ is more 
convenient than |q| and as usual u = [iq A — fi). All integrals are defined by Eq.(12). 

After the simple algebra has been performed within Eq.(56) our result is the equation of 
the fourth degree with respect to u;(£) 



2u 6 £ I B [ £ - 6(0 ] + u £ [ Ir (1 — b(0 2 ) -m 2 R - 2I K 



+ 2u £ I B [ I K ( 1 + d(0 ) 6(0 + i d(0 ( 2m il m - I K ) 
-/|£ 4 (l + rf(0 f + e d(0 2 ( 4 - 4m 2 Jl ) = 



(58) 



where = m(l — 31 z) is the renormalized fermion mass, and the functions and 6(0 
are given by 



rf(0 = | In ^4 

&(o=e-^(i-ain|^j 



(59) 



The dispersion relation (58) being very complicated is not solved analytically. However in 
the long wavelength limit (when £ — > oo) it can be simplified as 



to 2 - uj(I b - V m R) - (Ik - 2i]mI B ) ■ to 2 - uj(I b + r)m R ) - (I K + 2rjmI B ) 
and one finds the rather simple solution 

1 







w(0) 



rjm R + I B 



± 



[r] m R + I B 



+ (I K + 2rjmI B ) 



(60) 



(61) 



which demonstrates four well-separated effective masses: two of them pertain to quasi- 
particle excitations and other two to quasi-holes. Here rj — ±1, and the parameters m and 
H are nonzero. 
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